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Abstract 



We present a new method to derive low-lying N-dimensional quantum 
wave functions by quadrature along a single trajectory. The N-dimensional 
Schroedinger equation is cast into a series of readily integrable first order or- 
dinary differential equations. Our approach resembles the familiar W.K.B. 
approximation in one dimension, but is designed to explore the classically 
forbidden region and has a much wider applicability than W.K.B. . The 
method also provides a perturbation series expansion and the Green's func- 
tions of the wave equation in N-dimension, all by quadratures along a single 
trajectory. A number of examples are given for illustration, including a sim- 
ple algorithm to evaluate the Stark effect in closed form to any finite order 
of the electric field. 



PACS: ll.lO.Ef, 03.65.Ge 
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1. INTRODUCTION 



In this paper, we begin with the generalization to arbitrary dimensions of a new method 
in which the low-lying quantum wave function of a particle in a potential +V can be 
obtained by quadratures along the trajectory of a corresponding classical mechanical problem 
with —V as its potential. Let q = (gi, q2, - ■ ■ , Qn) describes the A^-dimensional coordinates of 
the particle, and $(q) is the ground-state (or low-lying) quantum wave function that satisfies 

i7<l>(q) = E^cl), (1.1) 

where 

H = -^V' + V{ci) (1.2) 
is the Hamiltonian of a non-relativistic particle of unit mass, and 

TV Q2 

V' = E#2- (1-3) 
We assume V{q) to be bounded from below; therefore we can choose 

V{ci) > 0. (1.4) 

Its minimum V{q) = may occur at more than one point, each with a non-vanishing curvature. 
As in [0, we introduce a scale factor g"^ by writing 

V{q)=gMq), (1.5) 

and consider the case of large g. We express 

$(q) = e-»^(^) (1.6) 

in terms of a formal expansion: 

^7S(q) = ^7So(q) + Si(q) + g-'S,{q) + g-'S^) + ■■■ (1.7) 



with a similar expansion for E, 

E = gEo + E, + g-'E2 + g-'E3 + ---. (1.8) 

Substituting ( |L6|) - ([T^ ) into the Schroedinger equation ([Tl] ) and equating the coefficients 
of (7^" on both sides, we derive 



(VSo)^ = 


2v, 






VSo • VSi = 


iv%- 






VSo • VS2 = 




- (VSi)^] 




VSo ■ VS3 = 




-2(VSi) 


■(VS2)]-^2, 



etc., where V is the A^-dimensional gradient vector whose components are V^ = d/dqi and 

i = l,2,---,N. (1.10) 

The underlying reason for considering the expressions ( |1.7|) and ( p..8|) is that for large g the 
ground-state quantum wave function has its amplitude centered around the minimum of V{q); 
near the minimum, V{q) may be approximated by a quadratic function of q, with the curvature 
of the corresponding potential surface proportional to g'^. Thus, for large g, the ground-state 
energy is 0(g) and so is the exponent of the ground-state wave function, as suggested by 
the ground-state energy and the corresponding wave function of a harmonic oscillator. The 
parameter g^^ serves as a measure of anharmonicity. 
It is useful to write the first line of ( p..9|) as 

^(VSo)^-^(q)=e, (1.11) 

where 



e = + . 
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(1.12) 



Equation ( |1.11|) is equivalent to the corresponding Hamilton-Jacobi equation in classical me- 



chanics, with —V{q) as the potential and e the total energy. To solve ( |1 . 1 1|) , we pick a 
minimum of V{q), say q = with 

V{0)=gMO) = 0. (1.13) 

Consider a trajectory q(t) which begins at q = when t = and ends at qT- when t = T; 
i.e., 

q(0) = and q(T) = qx- (1.14) 

Hamilton's action-integral is given by the minimum of the path integral 

T 1 ^ . 

So(qT, e) = y [- E - i-viqitmdt + T e, (1.15) 







i=l 



where 



(1.16) 



is the total time of the trajectory from the initial point to the final point qx- Since the 
potential in the classical problem = — f (q) < 0, the integrand in ( p..l5| ) is positive everywhere. 
Hence, given any final point qx and any total energy e > 0, the minimum of ( p.. 151) exists, so is 



therefore the classical trajectory. (For clarity, assume the minimum to be unique, a restriction 
that will be removed later.) As e — > 0+, the initial velocity q is proportional to e^/^ which 
also approaches zero; the total time T of the classical trajectory is oc e^^/^ oo, yielding 

lim T e = 0. (1.17) 

e=0+ 



The same limit e = 0+ of (pHSj ) gives So of ( plj) - ( CT ) 



To go back to the first line of ( p..9|) , we designate the end-point qx of the classical trajectory 
simply as q, the point of interest for the quantum solution $(q). Setting 



qx = q, (1.18) 
the solution So(q) of the first line in ( |1.9|) can now be written as 

T 1 ^ .2 

So(q)= J[^J:^^ -i-vmdt, (1.19) 



where the integral is along the trajectory qi{t) that satisfies the classical equations of motion 

q,= dv/dqi, (1.20) 

the energy conservation 

1 N 

^dtf ~ v{^{t)) = + (1.21) 



2 . , 

1=1 



and the initial and final configurations ( |1.14|) and (|1.18| ) 



The solution So(q) thus obtained is positive and continuous everywhere; it is an increasing 
function along the direction of the classical trajectory. In the case that the classical potential 
— t>(q) has other absolute maxima besides 0, say 

- t;(ai) = -v{a2) = ■■■ = -i;(a„) = 0, (1.22) 

So(q) is not analytic at these isolated points 

q = ai,a2, (1-23) 

For a total energy e > 0, along the classical trajectory the velocity must continue in its 
direction when the trajectory passes through any of the other maxima of — f (q), say q = an. 
Thus, as e 0+, VSq becomes zero at q = an, but its sign along the trajectory-direction 
remains the same before or after an; this forces VSq to develop a kink at q = an. At the 
initial point q = 0, while —v(0) = and therefore VSq = as e = 0+, So(q) is analytic, 
since different trajectories emanating from q = have to go along different directions. At 
infinity, it is easy to see that So(q) = oo, and therefore $o = e~^^" is zero. 



To solve the second equation in ([0|) , we require Si(q) to be also analytic at q = 0; this 
yields 

^0 = ^(V'So)at q=0. (1.24) 

It is convenient to consider the surface 

So(q) = constant; (1-25) 

its normal is along the corresponding classical trajectory passing through q. Characterize each 
classical trajectory by the So-value along the trajectory and a set of — 1 angular variables 

a = (ai(q), a2(q), • • • , aiv-i(q)), (1.26) 

so that each a determines one classical trajectory with 

Vaj ■ VSo = 0, (1.27) 

where 

j = l,2,---,iV-l. (1.28) 

(As an example, we note that as q — > 0, f(q) ^Y^^iQi snd therefore So ^J2^i<li- 

i i 

Consider the ellipsoidal surface So = constant = For / sufficiently small, each classical 
trajectory is normal to this ellipsoidal surface. A convenient choice of a could be simply any 
— 1 orthogonal parametric coordinates on the surface.) Each a designates one classical 
trajectory, and vice versa. Every (So, a) is mapped into a unique set (gi, g2, (In) with 
So > by construction. Depending on the problem, the converse may or may not be true; 
i.e. (gi, g2, Qn) (So, a) could be either one to one, or one to many. In the latter 

case, we regard the points in the q-space as specified by the coordinates (So, a), rather than 
the original coordinates (gi, g2, In)- 
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Write 

Si(q) = Si(So,«), (1.29) 

the second line of ( [1.9|) becomes 

(VSo)^(^)« = -V%-i?o, (1.30) 

and therefore 

Si(q) = S,(So, a) = / ^^[^V^So - Eo], (1.31) 

where the integration is taken along the classical trajectory of constant a. Likewise, the third, 
fourth and other lines of ( pT^ ) lead to 

Ei = ^[V2Si-(VSi)2],tq=o, (1.32) 
S2(q) = S.(So, a)=l ^^{^[V^Si - (VSOl - E,}, (1.33) 

^2 = ^[V2S2 - 2(VSi) ■ (VS2)]at q=0, (1.34) 

S3(q) = 83(80, a)=J (^{^[V^S^ - 2(V80 ■ (V82)] - E,}, (1.35) 
etc. These solutions give the convenient normalization convention at q = 0, 

8(0) = 

and (1.36) 
$(0) = 1. 

The above expressions for Eq, Ei, E2, ■ ■ ■ can be derived more directly by introducing 

f^(q) = ^8(q) 8o(q) 

= 8i(q) + g-'S^icO + g-'Ssiq) + ■■■. (1.37) 
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Thus, the ground-state wave function $(q) is 

$ = e-^^''-"", (1.38) 
and the corresponding Schroedinger equation ( |1.1|) gives 

gVSo ■ Va = ^[gV^So + VV - (V(t)2] - E. (1.39) 

At q = 0, f (q) = and VSq = 0, on account of (VSq)^ = 2v. Assuming Va to be regular, 
we have 

E = ^[gW^So + VV - (Va)2]3t (1.40) 

which leads to ([r24| ), (|L3^ ), (1.34), ■ ■ ■ for Eq, Ei, E2, ■ ■ -. 
Example. 

As an example, consider the special case 

nq) = y(?? + ?2 + --- + 4)- (1-41) 
From ( |1.9|) , one can readily show that the first line gives 

So(q) = |(g? + g2' + --- + g^) (1.42) 

and (|r^ leads to 

Eo = N/2; (1.43) 

the other equations ([T3I| ) - (1.35), etc., yield Ei = E2 = ■ ■ ■ = and Si = S2 = ■ ■ ■ = 0. 
The result is the well known exact answer 

$(q) = exp[-^{qf + qj + ■ ■ ■ + qj,)] (1.44) 

and 



E = ^N. (1.45) 

The organization of the paper is as follows: In Section 2, we discuss the case of a separable 
potential, and show that the solution of the wave function is a product of one-dimensional 
functions, as expected. The particular one-dimensional example of 

V{x) = ^{x^ - a^f (1.46) 

has been analyzed in ref. using our method. In this problem, besides the anharmonicity 
parameter g^^, there is also a much smaller barrier penetration parameter 

e-t5«'. (1.47) 

(The same problem has also been extensively studied in the literature[2-9], using mostly Feyn- 
man's path integration method. Quite often, the two parameters g and a in ( |1.46| ) are related 
by 2ag = 1, which makes V = |l/^(l — guY with y = x + a. The dimensionless anharmonicity 
parameter 1/ga^ then becomes Sg"^ and the barrier penetration parameter (|1.47| ) becomes 
In our (yf^^-expansion, we keep the parameter a fixed.) 
The ground-state wave function ipeven{x) is even in x, and the first excited state is an 
odd function ^oddi^)- Let £'e„en and Eodd be their corresponding eigenvalues. We can expand 
both their sum and difference as the following double series: 

hE,,, + E,,en) = E C^n(^7«')-"^e-t"^'^' (1.48) 

m,n 

and 

l-{E,,,-E,,,^)=4{-g'a')-2e-l^^' 

Z 71 

■Ecn.n(^7a3)-'"e-t("-i)^'^' (1.49) 

m,n 

where Cmn and Cmn are numerical coefficients, which can be explicitly expressed as definite 
integrals with the results 
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CmO = for all m, 

Coo = coi = 1, Cio = (1.50) 

„ _ 9 ^ _ 89 

<-2o - <-3o - etc. 

The terms with n = 1 represent the one-instanton contributions[2-9], n = 2 the two-instanton 
contributions, etc. As one can see from the derivation, an important tool, in addition to the 
expansion ( |L7|) - ( p..9|) discussed above, is the use of the Sturm-Liouville type of Green's 
functions. In this paper, we will show how to generalize such Green's function technique to an 
N-dimensional problem. In Section 3, we introduce a compact notation which will pave the 
way towards the generalization. 

Section 4 contains the main body of this paper. We depart from the initial approach based 
on the Hamilton-Jacobi equation; instead, we consider two N-dimensional Hamiltonians: 

iJ = _lv2 + y(q) (1.51) 

with, as before, 

V{q) = gM^) > 0, 

and 

n = H + eU{q). (1.52) 

Assuming that the ground-state wave function $(q) = e"^^'^''^ of the unperturbed Hamiltonian 
H is known, we set out to derive the perturbation series expansion for the corresponding 
ground-state wave function ^(q) of H with ell as the perturbation. Again, as we shall see, to 
each order in e, the perturbed \&(q) can be obtained by quadrature along a single trajectory 
connecting 0, the point that the unperturbed $(q) is maximum (i.e. S is minimum), to the 
point q of interest; the trajectory is determined by being the normal to the S = constant 
surfaces along the entire course of the trajectory. (It may be called the quantum trajectory, in 
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contrast to the classical trajectory determined by the Hamilton-Jacobi equation.) The form of 
the new perturbation series differs from the usual expression; the result is of course the same, 
as will be illustrated by a few examples. The new perturbation series formula will also be 
shown to cover the previous analysis in Sections 1-3 based on the classical trajectory derived 
from the Hamilton-Jacobi equation as special cases. 

In the same Section 4, the Green's function G of the N-dimensional Hamiltonian H will 
be derived, with G satisfying 

{H-E)G^\ (1.53) 

and E being the ground-state energy of H. Usually, the Green's function and the perturbation 
series expansion to each order e", require either an infinite sum over all excited levels of H, or 
equivalently, a sum over all possible paths through Feynman's path integration method. The 
Green's function and the perturbation series formula derived here are different; they depend 
only on quadratures along a single trajectory. 

In Section 5, we give examples to illustrate the applications of the new perturbation series, 
as well as to clarify the properties of the Green's function. In Section 6, we show how the 
trajectory-quadrature formulation can be extended to derive excited states, and in Section 7 we 
discuss how generalization can be made to allow for singular potentials, such as an attractive 
Coulomb potential with perturbations. As we shall see, with our new approach we can evaluate 
, e.g., the Stark effect to any finite order of the perturbation by quadratures along the radial 
trajectory. 

Although our methods are elementary, they become surprisingly powerful for the following 
reason. For example, in (1.38) by writing the wave function in the form of e~^^°~'^, where g is 
large, we cause the cross-term g'VSo • Vcr to be much larger than |[V^(7 — (Vcr)^]. Therefore 
the differential equation for a can be written order by order in such a way that |[V^(T — (Vcr)^] 
is known from previous orders, and VSq • Vcr = (VSq)^^ is inferred from the equation. In 
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effect, to each order, a now satisfies a first order ordinary differential equation (along the 
trajectory of constant a), not a second order partial differential equation. This makes possible 
a direct solution through quadratures. With variations, this theme runs through all the sections 
to be presented in this paper (sometimes e^"^ is written e^^, or and sometimes the known 
exponent is called gS, rather than gSo). 



For clarity, we consider the special case of dimension N = 2. The discussions below for a 
separable potential can be trivially extended to any N > 2. Let the two dimensional coordinate 
q in the previous section be designated 



2. SEPARABLE POTENTIAL 



q = {x,y) 



(2.1) 



and the quantum mechanical potential ^(q) and 



in (ITI) be 



V{q) = g^viq) = g^[v^{x) + Vy{y)] 



(2.2) 



and 



(2.3) 



As in (|1.4| ) and ( |1.13|) , we assume 



Vx{.x) > 0, 



(2.4) 



and at the point x = ?/ = 0, 



vM = ^y(O) = 0. 



(2.5) 



The classical Lagrangian density for the corresponding potential — f (q) is 
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l[x{tf+y{tf]+v,{x)+vy{y), (2.6) 
where the dots denote time-derivatives, as before. The equations of motion are 

X (t) = -tv^{x) 

and (2.7) 

y it) = i<y)- 

As in ( p..l4| ), we consider a trajectory that begins at x(0) = y(0) = (the maximum, or one 
of the maxima, of —v) at t = and ends at 

x{T) = XT and y{T) = vt (2.8) 

when t = T > 0, in which, xt and yT are variables independent of T. The action integral 
along the trajectory is given by the sum of 

A^{xt,T)=][\ x^ -{-v^{x))]dt 



and (2.9) 

Ay{yT,T) = ][\ y -{-vy{y))]dt 



The minimum Ax{xt,T) and Ay{yT,T), keeping xt, yr and T fixed, gives the Lagrangian 
equations of motion (|2.7| ), and the derivatives of A^ and Ay with respect to the final variables 

Xt, yr and T are 

£-A,{xT,T) =x (T) 
^AyivT^T) =y (T) 



|.A,(xT,T) = -e, (2.10) 



and 



■§pAy{yT,T) = -ey, 
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where and e„ are the energies in the x and y directions, given by 



1 -2 



and (2.11) 

. 2 

ey = ^y -Vyiy). 

For a given final configuration characterized by xt, Vt and T, both e^, Cy, and therefore also 
their sum 

e = e^ + ey (2.12) 

are determined. Following ( p..l5| ), we transform the independent variables from xt, Vt and T 
to Xt, Vt and e. We write 

So{xt, Vt, e) = A^{xt, T) + Ay{yT, T) +T e, (2.13) 

where T is a dependent variable, determined by 

T = T{xT,yT,e) = i^).,,,y,. (2.14) 

Given xt, Vt and e, the classical trajectory is given by the minimum of ( [2.13|) , from which T, 
Cx, Cy are also determined. Thus, we can regard 

ex = ex{xT,yT,e), 

and (2.15) 
Cy = ey{xT,yT,e). 



Hamilton's action integral ( p..l5|) becomes 



So{xT,yT,e) = Xo{xT,ex) + Yo{yT,ey), (2.16) 

where 
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Xo{xt, Ca:) = A^{xt, T) + T 6^ 



and 



(2.17) 



YoiVT, ey) = Ay{yT, T)+T a 



yi 



in which T is given by ( |2.14|) . Combining ( |2.10|) and (|2.17|) , we can readily show that 



f dSp \ _ (dYo\ _ (9Ay_\ _„ /rp^^ 



and 



(2.18) 



T. 



Along each classical trajectory, when the final time variable changes from T to T + dT, 
the end point moves from {xtiVt) to (xt+ x {T)dT, yT+ U {T)dT); correspondingly, 



dSo = (^)e. X {T)dT + (^)e, y {T)dT + e dT 



dxT dyx 
Combining the first two lines of (|2.18| ) with ( |2.19| ) and taking the limit e = 0+, we find 



(2.19) 



= {VSofdT. 



(2.20) 



Along the trajectory we also have 



dT = dxrl X (T) = dyr/ V (T), 



(2.21) 



which, together with ( |2.18|) , leads to 



dT = dxxl ( 



dyr/i 



dYn 



' dxT dyT ^ 

Substituting ( [2.20| ) and ( p.22|) into (|1.31J ) and labeling x-r, yr simply as x, y, we obtain 



(2.22) 



S,ix,y)=X,ix) + Y^iy), 



(2.23) 
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where 



Xi(x) = J dx{ 



dx ' ^2 (9x2 



with 



and 



and 



7? — l('^!2^o^ 

-^0,x 2\ (9x2 jat a;=0 



771 — l{d^\ 
-^0,y ~ 2 V (9y2 J at y=0 



(2.24) 



(2.25) 



-£^0 = + -£^0. 



(2.26) 



consistent with ( p..24|) . 

Likewise, S2, S3, • ■ ■ can be written as sums of the form X2{x) + Y2{y), X^{x) + Y^(y), 
■ ■ -, with 



v( ^ [, fdY, , I d^Y, dYr 



X,ix) = Jdxi^) {-^[^-^^^]-E.,h 



v( ^ f.l fdYo , I d'Y^ dY,dY2 



(2.27) 



where 



at x=Oy 
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El,y = 


1 
2 






E2,x = 


1 
2 




^ o o Jat x=0 

ox ox 


E2,y = 


1 
2 




^dY^dY^ 

^ o r^ Jat j/=0; 

ay ay 



etc.. The corresponding wave function factors e e e ^ ■■■ are all products 

^-g(Xo+Yo)^ g-(Xi+yi)^ e-<?-M^2+^2)^ etc.. 



3. A COMPACT NOTATION 

In this section we shall go over completely to the variables So, a and express the g^^ 
expansion (( |1.7|) and ( |1.31|) -(1.35)) described in Section 1 as iteration of an integral operator 
along a classical trajectory, that is for variable Sq and fixed a. As we shall see, our results can 
be written in a compact form by introducing an implicit matrix notation in which the indices 
are two values of S. 

In terms of the function (j(q), introduced by (|1.37| ), the wave function $(q) is 

$(q) = e-9So(q)-^(q)^ (31) 

Because 

(VSo)' = 2v, (3.2) 
the Schroedinger equation (pTlp is equivalent to the following equation for e~'^: 

[(^(VSo) ■ V + r]e-'^('i) = [-|v'So + E]e-''^^\ (3.3) 

where 

r=-^V^. (3.4) 
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Its solution is given by (|1.37| ) and ( p.. 401 ), which can be written in a more compact form, as 
we shall see. 

Consider the coordinate transformation 



Ql, 12, q3,---,(lN ^ So, «!, 02, • • • , «Af-l 



(3.5) 



with a = («!, 02, • • • , ctAT-i) denoting the set of — 1 orthogonal angular coordinates intro- 
duced in (g) - (pT^) . Define the Munction: 



^o(So — So) = < 



1 if < So < So 

if < So < So 



(3.6) 



where So and So vary from to oo. For So > 0, its derivative gives Dirac's (5-function: 

d 



dSo 



^o(So - So) = <5(So - Sc 



(3.7) 



Let us define Co{a) as a square matrix in the So-space, with its matrix element 



(So|Co(«)|So) = 9-%{So - So)/(VSo)2 



(3.8) 



An important feature is that each matrix element of Co{a) connects two points q(So,a;) and 
q(So,«) along the same classical trajectory that satisfies ( |1.20|) - ( [I.21J ), and therefore having 
the same angular variables ai, a2, ■ ■ ■ , cun-i- Throughout the paper, the variables Sq and Sq 
vary between and oo, with the end points and oo treated as boundaries. All gradients 
refer to differentiations with respect to the Cartesian coordinates; i.e., 

d 



and 



TV Q2 



(3.9) 



as before. On account of ( p.. 27]) and (|3.6|), for Sq > 0, 



VSo • V^o(So - So) = (VSo)'5(So - So), 



(3.10) 
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where the component of V6'o(So — So) is 

- d - 
Vj6'o(So — So) = 7^6'o(So — So) 

with Qj = qj{So,a). From (|3^) , we also have 

gVSo ■ V(So|Co(a)|So) = 5(So - So). (3.11) 

Like Co{a), Oq also denotes a square matrix whose matrix elements are 

(Sol^olSo) =^o(So-So). (3.12) 

Matrix multiplications follow the usual rule; e.g., 

/"OO — — — 

(SolCo^olSo) = / t/So(So|Co|So)(So|^^o|So). (3.13) 
Jo 

Introduce h^{Q) to be a diagonal matrix whose elements are 

(So|/ig(«)|So) = 5(So - So)/(VSo)^ (3.14) 

Thus, 

(So|^^o/io(«)|So) = ^^o(So - So)/(VSo)^ (3.15) 

but 

(So|/i^(a)0o|So) = ^^o(So - So)/(VSo)'. (3.16) 
Equations (|3.8| ) and( |3.11|) can be written in their matrix forms 

C, = g-%hl (3.17) 

and 

(7VSo-VCo = l, (3.18) 
20 



in which the a-dependences of Cq and are suppressed. 

Theorem 1. The solution of the wave equation (|3.3| ) - (|3.4| ) for e^'^ satisfies 

/•CO 

g-.(So,a)^l + / rfSo(So|[l + Co(a)T]-iCo(a)|So) 
Jo 

[--V2So + E]e-"(^0' (3.19) 

Proof (i) We first establish the existence of the inverse matrix [1 + Co{a)T]^^ . Assuming the 
opposite, there would be a column matrix f(a) which satisfies 

[l + Co(a)T]/(a)=0. (3.20) 



Applying VSq ■ V from the left, using ( [3.18|) and on account of T = — |V^, we have for the 
matrix element /(Sq, a) of f{a): 

gVSo ■ V/(So, a) - ^VV(So, «) = 0. (3.21) 

Hence, 

V ■ [e-25^«V/(So, a)] = 0. (3.22) 

Multiplying this equation by /(So, a) and integrating over all space, we derive, through partial 
integration, 

J ci^qe-2^^°[V/(So, a)]' = 0, (3.23) 

provided 

g~29Sojyj_Q 3^ So = 00. (3.24) 

It is clear that ( p.23|) implies V/(So, a) = 0; i.e., /(Sq, a) = constant, which makes 
-iVV = 0. It follows then from (ICTj ), 

/(So, a) = 0. 
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(ii) From (|3j), 



[giVSo) ■ V + T](e-'^ - 1) = [-^V'So + E]e-''. 



Next, applying ^((VSo) ■ V onto 1 + CqT and using ( |3.18|) , we find 



(3.25) 



(^( VSo) ■ V) [1 + CoT] = (7( VSo) ■ V + T; 



(3.26) 



therefore, ( [3.25| ) becomes 



[^(VSo) • V][l + CoT](e-'^ - 1) = [-|v2So + i?]e-^ (3.27) 



On the other hand, assuming ( [3.19|) we have 



e"- - 1 = (1 + Cor)-iCo[-|v'So + E]e-'' 



(3.28) 



which satisfies ( p.27|) . In deriving the above, we see that for / = e~'^ — 1, ( p.24|) is satisfied, 
which ensures the applicability of (1 + CqT)^^. The theorem is then proved. 

According to (|3.17| ), Cq = 0{g^^), and therefore, for g'^ >> 1, we may expand 



(1 + CoT)-'Co = Co + i-CoT)Co + {-CoT)^Co + 



(3.29) 



Neglecting 0{g'^), a = Si and approximating (1 + CoTy^Co by Cq, reduces to 



(3.30) 



Taking its derivative through gVSo ■ V, we have, as expected. 



^7[VSo-VSi]e-^^ = [VSo-i?]e-'^ 



(3.31) 



which, upon the approximation E = gEo, gives the second line of ([L9| ). Likewise, we can 
derive ( lOlD - (1.35) for Si, ^i, S2, ^2, ■ ■ ■ from (pl). 
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4. GREEN'S FUNCTIONS AND PERTURBATION SERIES 



As we have seen, the g^^ expansion enables us to probe the ground-state wave function 
$(q) along a classical trajectory; to each order g^", the A^-dimensional quantum wave function 
can be calculated through quadratures, by performing a finite number of definite integrals 
along the same trajectory. To the same order of accuracy, the corresponding energy is also 
determined by a finite number of differentiations at a single point where the potential has an 
absolute minimum, say at q = 0. This result may seem unfamiliar, since the Schroedinger 
wave function $(q) and its eigenvalue E depend sensitively on the boundary conditions at oo. 
But in our approach, the determination of $(q) is based only on a single trajectory connecting 
and q. The boundary condition at infinity appears almost automatically, because the classical 
action (i.e., Hamilton's action integral) satisfies the property: 

So(q) ^ oo as q ^ oo. (4.1) 

Since the classical action is calculated for a potential —V{q), which has its maximum at 0, this 
condition ( [^.1|) imposes only a general restriction on the potential; in any case, it is considered 
to be an input in the new method. To explore further the essential features of the underlying 
mechanism, we shall in this section generalize our formalism without any explicit use of the 
Hamilton-Jacobi equation. 

In this section we shall derive several new forms of the Green's function for the N- 
dimensional Schroedinger equation; these are all based on quadratures along a single specified 
trajectory. With the Green's function we will be able to arrive at a new perturbation series 
expansion, different from the usual familiar ones. As we shall see (near the end of the section), 
the material discussed in previous sections 1-3 will appear as special cases. As mentioned in 
the Introduction, the new Green's function will also enable us to extend the multi-instanton 
calculations developed for the one-dimensional problem in ref. [|l]| to arbitrary dimensions. In 
what follows, we shall first give a number of definitions, then for clarity, arrange the main 
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results in the form of three theorems (Theorem 2, 3 and 4 below). Examples of the new 
perturbation series and the Green's function will be given in Section 5. 
Consider two A^-dimensional Hamiltonians: 

H = -U/^ + V(q) (4.2) 

with V(q) > 0, as before, and 

n = H + eU{ci). (4.3) 

Their ground-state wave functions are $(q) and \l/(q) respectively: 

iJ<l>(q) = ^^(q) (4.4) 

and 

7^^(q) = ^^(q). (4.5) 

Let 

e^ = 8-E. (4.6) 

Write, as before, 

nq)=/^(q), (4.7) 

$(q) = e-^^^'i) (4.8) 

and 

^(q) = e-^'^('i)-"('i). (4.9) 

We have, from (0) - (p|), 
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2 

|-(VS)'-|v'S-/t; + E = (4.10) 

and 

^VS ■ Vr + ^[(Vr)2 - VV] = e{U - A). (4.11) 

In the subsequent discussions in this section, we assume that S(q) is known, and our objective 
is to compute T(q), assuming that 

|e| « 1. (4.12) 

We assume further that S(q) has an overall behavior similar to the Hamilton-Jacobi action 
So(q) given by ( |1.19| ); in particular, S(q) has an absolute minimum at q = 0. (This convention 
will be adopted here on.) Therefore 

VS = at q = 0. (4.13) 

(Note that from here on the absolute minimum of V{q) may or may not be at the same 
point.) Throughout this section, the role of the Hamilton-Jacobi action integral So in previous 
sections will be replaces by S. As in ( |1.26|) - (|1.28| ), we introduce a set of — 1 angular 



variables 



/? = (A(q),/?2(q),---,/?jv-i(q)), (4.14) 



which satisfy 



Vpj ■ VS = (4.15) 



with 



j = l,2,---,iV-l. (4.16) 



Each point q in the A^-dimensional space will now be designated by 
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(S,/3i,/?2,---,/?^-i), (4.17) 

AAA A 

instead of (gi, 52, ■ ■ ■ , ?Af)- At q, let S, /^a, ■ ■ ■ , form a set of orthonormal 
unit vectors, which are the normals of the S = constant, /3i = constant, ■ ■ -, Pn-i = constant 
surfaces at that point. Correspondingly, a line element d q can be written as 

- A ^1 A 

d q=S hsdS + J2 13 j hjd(3j] (4.18) 
i=i 

the gradient is given by 

A 



Aid ^"^ A 1 9 

V =S + E ^5-, (4.19) 



and 

T = -^V^ (4.20) 

can be decomposed into two parts: 

T = Ts + Tp, (4.21) 



with 



in which 



N-1 

h=I[ hj, (4.23) 



and 



It follows from (|4a8| ), 

/^i = [(vs)^ \ = [(v/5or ^ ■ ■ ■ , = mor ^ ■ ■ • • (4.25) 
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The volume element in the q-space is 

c/^q = hshpdSdp (4.26) 

where 

7V-1 

df3=X{ dpj. (4.27) 

The trajectory of a given /3 defines a continuous curve normal to the S = constant surfaces; 
along any such trajectory we will introduce the following matrices, 9, C{P), D{P) and G{P), 
each of which has its matrix element connecting two points (S,/3) and (S,/3) along the same 
trajectory. The matrix element of 6 is given by 

1 if < S < S 

if < S < S, 

We now define C, D and G in matrix notations by (cf. ( [3.13| )-( P7n| )) 



{S\9\S) =^(S-S) 



(4.28) 



C^g-V^, (4.29) 

D = -2e-^^ee^3^—ee-3^hsh0 (4.30) 
hp 

and 

G = D{l^TpD)-^ = {l + DTp)-^D, (4.31) 

in which the /3-dependence of the matrices is not exhibited explicitly and e^^^, hs, hp are all 

diagonal matrices. (The conditions for the existence of inverse matrices, like [l + TjsD)^^ will 
be given later, after ( [^.62| ).) For example, the matrix element of C is 

{S\C\S) = g-^e{S - S)/(VS)l (4.32) 



Similar to ( [3.13|) and ( p.l5| )-( p7l6D , the product of any two matrices, say A and B, defined 



along the same trajectory (i.e., sharing the same P) has its matrix element given by the usual 
rule 
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fOO — — — 

(SU5|S) = / dS(S|A|S)(S|£|S). (4.33) 

JQ 



The matrix element of its derivative -§^^3 is given by 

d . = d 



(S|^^S|S) = /~dS— (S|^|S)(S|S|S), (4.34) 
with ^ operating only on the left index S, whereas that of is 

(S|^AB|S) = (S|(^)B + A(-)|S). (4.35) 

Theorem 2. 

1. ^VS-VC = I, (4.36) 
2. D and G satisfy 

D = e-3^{l + CTsY^Ce^^ = e-3^C{l + TsCy^e^^ (4.37) 

and 

G = e-^S(l + CTy^Ce^^ = e-^SC(l + rC)-ie^s. (4.38) 
3. Furthermore, D and G are the Green's functions Ts-\-V — E and T + V — E: 

{Ts + V -E)D = \, (4.39) 

(i/-^;)^^ (T + y -£;)(? = I, (4.40) 

where I denotes the unit matrix whose matrix element is 5{S — S). 
Proof 

1. For S > 0, 

^^(S-S)=5(S-S). (4.41) 
On account of the orthogonality VS • VPj = 0, we have 

VS • V = (VS)2 A. (4.42) 
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Eq.flOeD follows. 
2. Define 



D = e^^De-^^ and G = e^^Ge-3^. (4.43) 
Then, from (4.31) we see that 

D = -2ee'3^^ee-^'^hshp, {AAA) 

^ = -2e'^^^ee~'^^hsh,, (4.45) 
oo tip 



and, because of (4.22), 



TsD = -p-e^9^ee-'''^hshp + 1. (4.47) 
Multiplying by C on the left, we find 

CTsD = 2e^e'^3SQ^-2gS^. +c = -D + C, 
hp 

and therefore 

(1 + CTs)D = C, (4.48) 

or 

D = {l + CTsy^C = C{l+TsC)-\ (4.49) 
Since according to (4.43) 

D = e-^^De^^, (4.50) 

(4.37) is proved. 
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By using (4.31), we derive 

dD 

as 



and 



d hpdD 



On the other hand, because $ = e^^^ satisfies (Ts + V — E)^ = and since 



we have 



g rr "9 ^hfs hp- 



2hshp dS hs hs 
Therefore ( |4.53|) leads to 



1 d hpdD. 



2hshf3 dS'hs dS 



+ {V-E)D = 1; 



i.e., D satisfies (|4.39| ) and is the Green's function of Ts + — -E. 
From ( |4.31j ) and the identity 



1 = [1 - TpD{l + TpD)-']{l + TpD) 



we have 



[Ts + V- E)D = [1 - TpD{l + Tf,D)-\l + TpD), 



I.e. 



{Ts + V - E)D{1 + TfsD)-' = 1 - TpD{l + TpDyK 
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From ( |OlD , D{l+TpD)-^ is G; we derive 

{T + V - E)G = 1, (4.60) 

which establishes (4.40). 

From G = e^^Ge^^^ and since e^^ commutes with Tg, it follows that 

G = D{1 + Tf}Dy\ (4.61) 

By using D = G{1 + TsG)-\ we have 

l + TpD = l + Tf,G{l + TsG)-' 

= {l + TsG + TpG){l + TsC)-\ 

Its inverse is 

(1 + T^D)-^ = (1 + TsC)(l + TC)-i 

which leads to 

G = D{1+ TsC)(l + TC)-^ = C{1 + TC)-\ (4.62) 
and that gives (4.38). From the above expression, we see that 

G~^ = G-^ + T. (4.38)' 
Likewise, from (4.49) and (4.61), we can write 

D'^ = G-^ + Ts (4.37)' 

and 

G'^ = D'^ + Tp. (4.31)' 
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These relations make transparent the connections between the original equations (4.38), (4.37) 
and (4.31). 

To complete the proof, we have to examine the conditions, under which 1 + CT , 1 + CTs 
and 1 + DTp have inverses. 

As in ( [3.20|) , we assume /i = /i(S,/5) to satisfy 



(i + cr)/i = o. 



(4.63) 



Operating on its left by —2gVS ■ V, we derive 



- 2gVS ■ V/i + Wi 



0. 



Hence 



V ■ (e-'ssyj^) _ 0. 



(4.64) 



Multiplying by fid 



q and integrating over all space, we have 




(4.65) 



provided 



e-'^'/iV/i = 



at 



CO. 



(4.66) 



From (|4.65| ), we see that the only solution is V/i = 0; i.e., /i = constant. However, from 
the original equation (|4.63| ), it follows then fi = 0. Consequently, (1 + CT)~^ exists provided 
( F^ ) holds. 



Next, we examine the equation 



(l + CTs)/2 = 0. 



(4.67) 



Operating gVS ■ V on its left, we find 
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(7VS-V/2 + Ts/2 = 0. (4.68) 
Multiplying (|4.68|) by f2hshpdS and integrating from S = to oo, we have 

^(^^-^"i"'^-''' («^) 



Assuming 



{f2^^)e-^'^hp = Q at S = and S = oo, (4.70) 



(which is the S-component equivalence of ( |4.66|) ), we derive 

Recognizing ci^q = hshpdSd(3 is positive, it follows then 8/2/88 = 0. Substituting this back 
to ( 14.671 ), we find 

/2 = 0. (4.72) 

Hence, (1 + CTs)-^ exists, provided ( F70| ) holds. 
Lastly, we consider 

(l + DT^)/3 = 0. (4.73) 
Applying Ts + V — E from the left and using (|4.39| ), we obtain 



{Ts + Tfs + V-E)fs = 0. (4.74) 

Thus, the inverse of (1 + DT^) exists, provided we restrict ourselves to a Hilbert space which 
excludes the ground-state of H. Theorem 2 is then proved. 

It is important to note that the existence of (1 + CT)"\ and therefore G, as defined by 
(4.38), is free from the last restriction. 
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We return now to (|4~^ )- (|4.6| ). The following theorem will enable us to derive the ground- 
state wave function \E' = e"^^""^ of H + eU in terms of e^^^, the ground-state wave function 
of if: 



Theorem 3 



e"" = 1 + (1 + CT)-^Ce{-U + A)e~\ (4.75) 



or equivalently, 



^ = e-^s + (1 + DTp)-^Dt{-U + A)^ 

= e^s^ + Ge(-[/ + A)^. (4.76) 

Proof . Eq.( |4.11| ) can also be written as 

(^VS ■ V + T)e-^ = e{-U + A)e-\ 

where T = -^V^, as before. Since from ( |06D ^VS-VC = 1, and since V(e-^-l) = Ve"^, 
the above expression gives 

^VS • V(l + CT){e-^ - 1) = e{-U + A)e-\ (4.77) 

Hence, 

e"^ - 1 = (1 + CT)-^Ce{-U + A)e-\ (4.78) 



as can be derived by substituting ( |4.78|) into i ^.71\ ). Thus ( |4.75|) follows. Multiplying (K.78| ) 
by e~^^ on the left gives 

^ = e-f^ + e-^^(l + CT)-^Ce{-U + A)e-^ (4.79) 



which leads to (|4.76| ), on account of (|4.31j ) and (4.38). Theorem 3 is thereby established. 



From (|4.75| ), we may expand e ^ as a power series of e: 
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e-" = 1 + e(l + CT)-^Ci-U + A) + + CT)-^C{-U + A)]^ 

+ • • ■ + e"[(l + CT)-^C{-U + A)]" + • • • . (4.80) 

Likewise, from (|4.76| ), we derive the perturbation series 

^ = ^0 + e^i + e2^2 + • • • + e^-^n + • • • , (4.81) 

where 

^0 = e"^^, (4.82) 

and for n>l, 

^„ = [G{-U + A)]"e"f^, (4.83) 

with G = (1 + DTpY^D, or e^^il + CT)-^Ce<^^. 
From (|4.83| ), we have 



^„ = G'(-f/ + A)^„_i. (4.84) 

Eqs. ( |4.76| ) and ( |4.83|) can be directly established by using {H - E)G = 1; i.e., G being the 
Green's function of (H — E). Likewise, since 

[Ts + V - E)e-^^ = 

and 

(Ts + V - E)^ = [-Tp + t{-U + A)]*, 

\1/ also satisfies 

^ = e-s^ + D(-T^ + e(-f/ + A)]^, (4.85) 

on account of ( [4-.39| ), D being the Green's function of (Ts + V — E). 
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The new perturbation series (|4.81| )- (|4.83| ) is most effective, when it is complemented by 
an additional expansion in powers of g^^, the anharmonicity parameter. From C = g^^Oh"^ 
and (4.37)-(4.38), we see that the three Green's functions C, D and G are all 0{g^^). The 
n*'^ order perturbative wave function e"\l/„ is given by (|4.83|) . By using (4.32) for G, \&„ can 
be written as a power series in D, and therefore also in g^^: 



^„ = {[D - DTf,D + (DTpYD + {-nDTp)"'D + ...](-[/ + A)}"e-^s^ (4 gg^ 



a double definite integral along the trajectory. Through (4.38), can also be expressed as a 
different power series in g^^ in terms of G: 



in which each C is a single definite integral along the same trajectory. The equivalence between 
these two expressions rests on (4.37), which gives 



The identity between this power series expression of D and its double integral form given 
above has its origin in the Sturm-Liouville construction of the Green's function for a one- 
dimensional Schroedinger equation. Recalling that e^^^ is a solution of the second order 
ordinary differential equation in S, 



in which each D assumes the form given by (4.31); i.e.. 




■qr^ = e-^^{[C - GTG + ■■■ + {-)'^{GT)'^G + ■•■](-[/ + A)}", 



(4.87) 



D = e-^^\G - GTsG + ■■■ + (-)™(CTs)'"C + ■ ■ -Jef^. 



{Ts + V- E)e-^^ = 0, 



(4.88) 



we can form an irregular solution F which satisfies the same equation 



{Ts + V- E)F = 



(4.89) 
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with 

s 



F(S) = e-^^ / e^^^(|^)3, s^S, (4.90) 

in which the angular variable (3 is kept fixed and not exhibited explicitly in the argument. The 
Green's function D is related to the Wronskian-type expression by 

(S|D|S) = 2[e-^^F(S) - F{S)e-^\hshp)^^pe{S-S). (4.91) 

In ref. (jH], we used a one-dimensional example to demonstrate how this form of the Green's 
function can enable us to evaluate the multi-instanton expansion. Here we have extended the 
application of a one-dimensional Green's function to an iV-dimensional problem. 
Theorem 4. The energy shift eA is given by 

where d^<\ = hsh^dSdp, as in ( [4.26| ) and the integration is extended over all space. 
Proof From (4.31) or (4.91), the matrix element of D{(3) is 

where the subscripts S,/5 and S,P indicate the arguments of the functions inside the paren- 
thesis. 

Let X denotes the column matrix 



iS\DmS) = -2e-^^ £ dS{^)=/{S - S) ■ e'^%-^^{hsh,)^^,, (4.93) 



x(S,/3) = [-T^ + e(-f/ + A)]M/(S,/5). 

Eq.( pl-.85|) becomes 

^ = e"^^ + Dx- 

As S — > oo, the left side should be zero, but the right side is, on account of (|4.93| ), controlled 
by the S-integration in the upper region when S is near S; i.e., writing Dx as the product of 
-2e"^^ times 
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we see that the factor to the right of (and including) the S— integration sign approaches e^^^ 
times 

/ dSe-^^{hshp)^^px{S,P), (4.94) 
which would make Dx oo unless ( H-.94| ) is zero. Thus, the convergence of \1/ at oo requires 

POO 

/ dSe-^^hsh3[-Ta + e{-U + A)]^ = (4.95) 
Jo 

at all p. Integrating over the angular variables d(3 = Ufji d(3j, and because of 

-h.h,T, = -^-^^^ (4.96) 

and 

' dphshpT^m = 



we derive (|4.92| ). This result is well known. (Since (iJ + ef/)\l> = (i? + eA)\l>, the multipli- 



cation of e"^^ on the left and integrating over all space gives immediately (|4.92| ).) We have 
demonstrated that the same result can also be derived, though somewhat awkwardly, by the 
Green's function method developed in this paper. 

The usual perturbation series requires to each order n, either a product of infinite sums 
over all excited levels of the unperturbed Hamiltonian H, or a sum over all possible paths 
through Feynman's path integration method. The perturbation series formula derived here 
is different; it depends only on quadratures along a single trajectory, as will be illustrated by 
examples in the next section. 
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Remarks. Before leaving this section, we wish to address the connection between the 
perturbation series developed in this section and the previous series expansion given by ( |3.28| )- 
( p.29| ). In sections 1-3, we discuss the relation between the Schroedinger wave function 



and the classical Hamilton's action integral Sq, which satisfies 

g\VSof = 2V = 2gh. 

Thus formally we may regard e~^^° as the solution of a Schroedinger-like wave equation whose 
eigenvalue happens to be 0: 

(_lv2 + \4)e-3^° =0 (4.97) 

where 

Vc = g'v - ^^^V^So. (4.98) 



Since e satisfies 



(-^V^ + V; + ^gV'So)e-^^ = Ee'^^ (4.99) 



we may treat (|4.97| ) as the unperturbed Schroedinger equation, and ( [4-.99|) as the perturbed 



one, with ^gV'^So as the perturbative Hamiltonian and E as the energy shift. The corre- 
sponding small e-parameter that characterize the perturbation becomes g^^, since ^gV'^So is 
0{g~^) times Vc- This explains the similarity between ( |3.28|) for e^"' and ( |4.75|) for e^^; it 



also identifies the g -expansion with the perturbative e-expansion. 
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5. EXAMPLES 



From ( |4.31| ) and (4.38), the Green's function G can be either expressed in terms of the 
single integral form C , or equivalently, the double integral form D: 

G = e-^^{l + CT)-^Ce^^ (5.1) 



or 



G={l + DTp)-^D. (5.2) 

These two complement each other in their applications. As in (|4.8| ) - ( |4~9|) , let e~^^ be the 
ground-state wave function of H . The corresponding ground-state wave function 

m = e-^^-^ (5.3) 

of if + ef/ can also be expressed either in terms of C (eq. (|4.75| )): 

e-^ = l + {l + CT)-^Ce{-U + /\)e-\ (5.4) 
or in terms of D (eq.( |4.76|) ): 

= e-^^ + {I + DTp)-^ Dt{-U + A)-^ . (5.5) 



As shown in ( |4.80|) and (|4.86| ), both can be used for the perturbation series expansion. In this 



section, we shall examine several examples which illustrate the different merits of these two 
formulations. 

For clarity, consider the one-dimensional case, in which ( |5.2|) and ( |5.5|) reduce to G = 

and 

m = e-^^ + De{-U + /\)^! . (5.6) 

As in (|4.43| ), introduce 
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D = 



which, in accordance with ( [^.48|) and since Tg = T in the one-dimensional case, is related to 
C by 

(1 + CT)D = C. (5.7) 

In terms of D, (|5^) becomes 

e-^ = 1 + De{-U + A)e'\ (5.8) 

In order that the above expression is equivalent to the alternative form ( |5.4|) in terms of 
C, we need to convert ( p7|) into 

D = {1 + CT)-^C. (5.9) 

According to (^^BSj), this requires 

e-'^^/iV/i = (5.10) 
at the boundary, where, for (|5.8| ) , the function fi is 

fi = D{-U + A)e-\ (5.11) 
and in one-dimension, the boundary consists of 

x = ±oo. (5.12) 



We recall further that C satisfies ( |4.36| ) which is now an ordinary first order differential equation 



dSdC 

ax ax 

this determines C up to an integration constant, which can be set by requiring 

^ = at S = 0. (5.14) 
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Now, D is the solution of the second order differential equation 

. 1 



,,^.+V'E)D^l, (5.15) 



and therefore contains two independent integration constants. In view of ( |5.9|) and (|5.14| ), we 
require 

D = at S = (5.16) 
and write for the right hand side of (|5.11j) 

X y 

D{-U + A)e-^ = -2 J e^^^^y^y J e~^<^^^'\-U + A)e'^^''^dz (5.17) 

-oo 

SO that ( |5.10D holds at the boundary x = — oo. (For convenience, we choose S = at a; = 0.) 
In order to satisfy ( |5.10| ) at the other boundary x = oo, we set 



oo 



Otherwise, ( [5.17|) would ~ e^^^ at a; = oo, violating (|5.10| ) - ( |5.11|) . Thus, the energy shift 



A is also determined, and ( |5.18| ) is a special case of ( |4.92|) . On the other hand, if we use 
( |5.1|) , expressing G in terms of the single integral form C, then A would be determined by 
differentiations at S = 0, similar to (|1.32| ) and (1.34). In the following, we shall solve two 
examples by using G, first in terms of C, then in terms of D. 
Example 1. Consider a one-dimensional harmonic oscillator with 

*-4^ + T^^ (5.19) 

and a perturbative potential 

eU = ex^P, (5.20) 

where p is a positive integer. The unperturbed ground-state wave function is 
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Since S = \x^, we have 



hi = (dS/dx)-^ = x-^ = {2S)-\ 



For n > and in accordance with ( |5.14| ), 



Jo g2S g{2n) 



Lemma 1 



m=l 

where 



2n-l 

i n-l 



g2n' 2g^{2n-2y 

(2n - l)(2n - 3) •• -(2171 + 1) , 
r = .1 i . L for 1 < m < n — 1 

mi2g)"+^ iSmSri 

and in particular 

(2n- 1)!! 



Proof Since T = -|^, 



(1 + CT)-'C = C- CTC + {-CTyC + ■ ■ ■ + (-CT)™C + 



and 

2g 



TCx^"" = —{2n - l)x2"-2, 



we find 



^^^"^ ~2g^{2n-2f " ^ ' 

/^T^2^„2n _ (2r2 - l)(2n - 3) 2n-4 _ p ™2n-4 

'^^^^"^ " 2V(2n-4) ' 
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• ■ -, until 



n~ln^2n ^ (2n-_l)!! ^2 _ -n ^2 



• 2 _ p 2 



But, (-CT)"Cx2" = CTi „, which by itself would be oo. On the other hand, 



Lemma 1 is then proved. 

As in ( f^.9|) , the eigenstate of 



1 '^^ 5'^ 2 2r, 



is e-^'J'^'-r Write 



e ^ = 1 + ^a^x^'. 

From ( |r75| ) 

oo oo 

azx^' = (1 + CT)-iCe(-a;2p + A)(l + a^a;^'). 
1=1 1=1 

By using ( ^.24|) and ( |5.27|) , we derive 

oo 

1=1 

and for n > 1 

oo 
1=1 



where n = for m > n; otherwise it is given by (|5.25| ). From ( |5.29|) ,we see 



eA = — oi- 



Expand 
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and 



eA = eA(l) + e2A(2) + e3A(3) + 



an = ea„(l) + e^an{2) + e^an(3) + 



it follows then 



A(l) = Ti p, a„(l) = -r„ p, 

p 

^(2) = X/ rn+p + Ti p) (5.32) 



m=l 
V 



a. 



i(^) ^ y m n m+p mJ^l p) ) 

m=l 



etc.. (Because of ( |03| ), a„(2) = for n > 2p.) 
As a special case, for 

ef/ = ex^ (5.33) 

and therefore p = 2, we find 



.A(l) = ^. .^A(2) = -^, (6.34) 



the same as the usual perturbation results. 

Example 2 Using the same one-dimensional harmonic oscillator Hamiltonian ( [5.19| ) as the 
unperturbed H, we consider now an odd power perturbation potential 

eU = ex^P^\ (5.35) 

The analysis given below is parallel to that in the first example, but with some changes, as 
will be indicated. 
For n > 0, 
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It is straightforward to establish the following: 
Lemma 2 



where 



(1 + CT)-'CX^'^^' =T.lrn nX^"^^' (5.37) 

m=0 



1 n 



Tn n /„ -| \ ) Tn— 1 



g{2n + iy g\2n-iy 
n{n — 1) ■ ■ ■ (m + 1) 

gn-m+1^2m+ 1) 



0<™<n (5.38) 



and in particular 

n! 

70 n = 

Instead of ( |5.27| ), we now set 

oo 

e"^ = 1 + ^ hnx"". (5.39) 

n=l 

From ( |475D , 



X] bnx" = (1 + CT)-iCe(-x2p+i + A)(l + J2 hnx""). (5.40) 

n=l n=l 

It is convenient to extend the definitions of 7^ „ and « (as in (|5.25| )), by defining 

7m n = n = for TU > u (5.41) 

and 

To „ = 0. 

Using the two Lemmas, we derive 

A = Xl(^2«+i Ti i+p+i - A b2i Ti i), 
I 

hn = -e^(62«+i r„ i+p+i - A b2i r„ ;) (5.42) 
I 
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and 

hn+i = -e7n p ~ eXl(^2« In l+p - A b2l+l 7n «)• 
I 

Here and in the following, whenever the range of the summation index is not exhibited, it 
reads 

/ (orm, n) =0, 1, 2, ■ ■ ■ , oo. (5.43) 
By comparing the first two equations in ( [5.42|) , we have 

eA = 62- (5.44) 
It can be readily verified that the following power series expansions hold: 

eA = e2A(2) + e^A(4) + t^^{<o) + ■■■ 

h2n = e%n{2) + e%2n{4:) + e%2n{Q) + " " " , (5.45) 
b2n+l = e&2n+i(l) + e^fe2n+i(3) + e^62n+i(5) H , 

in which 

&2n+l(l) = -In p, A(2) = -^7; pTi l+p+l, 

I 

^2n 

(2) = E 

ll p^n l+p+li 

I 

&2n+l(3) = - ^ll p^m l+p+lln m+p - A(2) ^ 7„ rnlm p, (5.46) 
l,m rn 

A(4) = -A(2)5:(7, 



Lm 



n. m+p7l n+p+1) 



etc. 



Z,m,n 

In the special case of p = 0, we have 

eU = ex, (5.47) 
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the above expressions lead to 



e^A(2) = -— (5.48) 

and all other A(2n) = 0. This confirms the exact result, since 

V + eU = ^-^x' + ex = ^{x + - ^. (5.49) 

The lowest eigenvalue of is |, and that of + ef/ is | — Likewise, ( |5.45| )- (|5746| ) yield 

hi = ebi{l) = -e7o o = --, 

9 
e2 

b2 = e%i2) = e%o o^i i = (5.50) 



h = e^&3(3) = -e^7o o7i i = ^ 



3 



6g^ 



etc., so that 



e"" = 1 + ^ 6„x" = e"^/^ 

n=l 

which leads to the expected expression 

e-'^-' oc e"^(-+^)'. (5.51) 

So far we have only used the Green's function G in terms of the single integral form C. 
Next, we will repeat the same examples, but using the double integral form D. 
Example 1 (repeat) We return to the same simple harmonic oscillator problem ( |5.19| ) - 
( F20| ), with 

I 9' 2 

H = — + —x^ 

2dx^ 2 

and the perturbation 

eU = ex^P. 
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In order to derive the perturbation series ( |5.27D - ( |5.32D for e ^ and A using D, we have to 
deal with integrals of the type 

X y 

W = -2 J e^y'dy J e-3''z^''dz. (5.52) 

-oo 

Introduce 
( |5.52| ) becomes 



Dx2" = — —Je"'dr] J e-^'C'"f^C- (5-53) 



-oo 

Let Hi{Q be the usual Hermite polynomials, obtained by the generating function 



oo J./ 

e-^+^*^ = EU(C). (5.54) 

1=0 



The integrand in ( |5.53|) can be expressed in terms of Hi{() through 



For / > 1, 



therefore 



r-t y„,„,gf:„„), g^.-.4C). (5.55) 



>3 • I 

J e^'df] J e-<'Hi{C)dC = -J Hi.Mdv 



iJKO-^KO)] (5.56) 



-oo 

1 



2V 



and 



2 ' 



DHii^x) = --J e^'dr^ J e-^' Hi{C)dC 



-oo 
1 



[HiiO - HiiO)]. (5.57) 
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Since -f^o(C) = 1. by separating out in (|5.55| ) the term m = n in the sum, we have 

^0 22"m!(2n - 2m) 



= „ + E ,^,, ^2.-2^(0, (5.58) 



where 

(2n-l)!! 
(2^7)^^ ' 

the same as (|CT|) . From (|537|) and (|53^ ), it follows that 



The identity between the above expression and (1 + CT)^^C(x^" — Fi „) given by ( |5.24| ) 
can be established by first applying T onto ( |5.57| ). Through the well known formula 

^,m) = 2^^^m)-2im), (5.60) 

we find 

TDHii^x) = ^[-x^Hii^x) + IHii^x)]. (5.61) 



Next, apply C onto (|5.61| ) 



CTDHii^x) = -j —TDHii^x) 





^m^x) - + CHii^x). (5.62) 



Combining (|5.57| ) and ( |5.62|) , we obtain 



(1 + CT)DHi{^x) = CHii^x). (5.63) 
Hence, the multiplication of (1 + CT) and ( |5.59|) leads to 



(1 + CT)D{x''^ - Ti „) = C{x'^ - Ti „), (5.64) 
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which gives the equality between (|5.59| ) and ( |5.24|) . Thus, both forms of G, (|5.1| ) and ( ^^ , 
yield the same perturbation series expansion ( |5.31| ) - ( |5.32|) . 
It is instructive to note that in (|5.24| ), 



n 



(1 + CTY^Cix^^ - Ti „] = ^ „a;2-, 



m=l 



each term, (1 + CT)"^Cx^" or (1 + CTy^CTi „, is by itself oo, which necessitates the 
combination that appears on the left-hand side. Using the single-integral form C for G, the 
energy shift A is required to make the integral finite at x = 0. On the other hand, Dx^'^, or 



X y 



DT, n = -2j e^y'dy J e-^^'li ^dz, 

-oo 

is well-defined at any finite x. As x oo, either Dx^" or 

oo 

DTi n ~ -2[ri „ J e~<^''dz]xe^^' (5.65) 

— OO 

1 2 

approaches oo, which, when multiplied by e~29^ , would lead to a perturbed wave function 
g-gs-r divergent at a; = oo. As shown in (|5.17| ), this requires the energy shift A to be 
determined by ( |5.18| ); in this example, it is equivalent to have the subtraction —DTi „ in 

im)- 

Example 2 (repeat). To complete the one-dimensional analysis, we continue with the same 
simple harmonic oscillator H as the unperturbed Hamiltonian, but with the perturbation given 
by 

eU = ex^P^K (5.66) 



Since by repeated partial integrations, 

J 2(7 o 



-co 

g2 gn 
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we find 

X y 

-oo 

n 

/ , Im ) 
m=0 

the same as ( |5.37| ) for (1 + CT)"^Cx^"+\ with 7^ „ given by (|5.38| ). It is straight forward 
to see that using either D, or (1 + CT)^^C , we can arrive at the same perturbation series 
(131) and ( p^ . 
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6. EXCITED STATES 



In this section, we apply the new approach to excited states. 
Let e"^^ be the ground-state of (1.1) 

where 

and, as in (1.4)-(1.5), 

v(q) > 0. (6.1) 

In this section, we assume e~^^ to be already known; our purpose is to derive the excited 
states ^ex, which satisfies 

H^ex^{E + S)^ex (6.2) 

Write 

= X(q)e-^^ (6.3) 

Since 

V'^ex = (V'x - 2^Vx • VS)e-^s + x(V'e-^S), 

(6.1)-(6.3) lead to 

^VS • Vx - l^'x = ^X- (6.4) 

The expansion (1.7), 

gS^gSo + Si + g-'S2 + g-^S3 + ---, (6.5) 
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will now be accompanied by similar expansions for x 3nd £: 

X = Xo + 9~^Xi + 9~^X2 H (6.6) 

and 

£^g£o + Si + g-'£2 + ■■■■ (6.7) 

Substituting (6.5) - (6.7) into (6.4) and equating the coefficients of g~'^ on both sides, we 
obtain the following first order partial differential equations for xo: Xi: X2, • • 

(VSo ■ V - £:o)xo = 0, (6.8) 
(VSo • V - £o)xi = (-VSi • V + ^ + £^)xo, (6.9) 
(VSo ■ V - So)x2 = (-VSi . V + + Si)xi + (-VS2 ■ V + S2)xo, (6.10) 

etc. To see how these equations can be solved we first give a simple example and then address 
the general solution. 
Example Take the example 

^(q) = li^hl + ^hl + • • • + ^Wn). (6.11) 



and correspondingly, 



So(q) = + ^2^2 H \- J^nQn) 

Si(q) = S2(q) = --- = 0, (6.12) 



|(i/i + i/2 + --- + i/^). (6.13) 



Introduce, as in (1.26)-(1.28), the variables 

(So, a) 
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with a — (cKi, 0:2, • • • , ccjv-i) and 

VSo • Vaj = 0. (6.14) 
In terms of these new coordinates So and a, (6.8) becomes 

(VSo)^(^). = E,, (6.15) 

which will be integrated along the constant-a trajectory. The constant-a trajectory is deter- 
mined by the solution of the Hamilton-Jacobi equation (1.11); it describes a classical trajectory 
with —V as the potential and a positive infinitesimal energy, starting from q = at time near 
—00, and q = (?i, ?2, • • • , Qn) — (So, 01) at time t. As t increases to t-\- dt, the end of the 
trajectory moves from q = (So, a) to 

{qi + dqi, q2 + dq2, ■■■,qN + dqN) = (So + dSo, a) (6.16) 

keeping a constant. In accordance with (2.20) - (2.22), we have 

ciSo dqi dq2 dqN . „ 

dt = = = = = . 6.17 

(VSo)^ uiqi ^2^2 J^nQn 

When So 0, each q^ must also —>■ along the trajectory. Since Xo is a single-valued function 
of q, we can classify xo by its power dependence on q^, as q^ 0. Take this to be 

Xo^qi^qT-'-lN"- (6.18) 

Since 

From (6.17), we see that along the classical (constant a ) trajectory 

= (VS^' ^'-'^^ 
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and consequently 



' dSo (VSo)2 ^ ' dlnq, 

As So — > 0, under the assumption (6.18), we have 



E-^^- (6-21) 



therefore (6.15) yields 



(VSo)^(^)«-En,^.; (6.22) 



N 

^0 = XI (6.23) 

1=1 



and 



(VSo)^(^)a = E^^^^ (6-24) 

along the entire trajectory for all So > 0. (It is reassuring that (6.23) is independent of a.) It 
follows then 

Xo = qrQr---Q7 (6-25) 

at all q. 

In this example, because Si = 0, (6.9) becomes 

(VSo)^(^)a - Soxi = ^iXo + ^V^xo, 
which, on account of (6.15) and (6.25), can also be written as 

Xo(VSo) [^(— )]a = ^1X0 + ^[ -2 + -2 + • • • + -2 ]Xo- (6.26) 

Thus, keeping a fixed, 

/■ C^Sq 1 711 (711-1) , 712(^2-1) , , nNjuN-l) ^^ , . 

^ y (vs^^^^ + 2^-^- + + • • • + — — ^'-'^^ 
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where the integration is along the (classical) trajectory normal to the So = constant surfaces. 
We leave (6.27) in the indefinite integral form, since the integration constant can be eliminated 
by the transformation 

Xi ^ Xi + constant • Xo (6.28) 

which affects only the overall normalization factor of x- (This convention will be adopted 
below for the general case as well.) 
For any £ 7^ 0, and any partition 

£1 = rriiVi + m2Z^2 H \- rriNi^N 

(where rrii can be an arbitrary number) gives 

/ j^,S, = m.lng, + m.lnq, + • • • + m^lng^. 
J (VSo)^ 

This leads to an inadmissible wave function, and consequently 

£1 = 0. (6.29) 
Hence, by using (6.17), (6.27) and (6.29), we find for this example 



Xo r ^iK - 1) , - 1) . . njv(njv - 1) 

^^""T^ — — + — — + + — — ^^-^^^ 



Recalling that the Hermite polynomial Hn{x) is 



the sum 



X = Xo + ^ ^Xi + ■ 



can be shown to be, apart from an overall normalization factor, 
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Hni{V9i^qi)Hn2{\^q2) ■ ■■Hnj.i^/gv^qN)- (6.31) 

Returning to the general case, because q = is the minimum of v(q), as q — > 0, i'(q) 
depends quadratically on q. Write 

^(q) ^ + ^ql + • • • + ^Wn)-. (6-32) 

correspondingly, 

So(q) + ^2ql + h i^Nq%)- (6.33) 

Again, we shall classify Xo(q) according to its behavior as q — > 0, by using (6.18). As in this 
example, this leads to 

So = nivi + n2i^2 H h uni^n, (6.34) 

where, as before, ni, ^2, ■ ■ ■ ■, are positive integers. From (6.8), it follows that 

^"^° = ^°/(^ ^^-^^^ 
with the integration taken along the trajectory of constant a, and the integration constant is 
determined by the normalization condition (6.18), as So 0. 

It is convenient to characterize the limit So 0, by introducing along the trajectory an 
overall scale factor A (of the dimension [^j]). Because of (6.17), for sufficiently small Sq, we 
set 

Sq oc A^, qi cx A, (6.36) 
and therefore, on account of (6.17), 

.^^°.„ = — dhiqi = — (ilng2 = • • • = — dinger oc dlnX. (6.37) 

(VSo)^ V2 UN 

We now turn to (6.9); as in (6.26)-(6.27), it can be written as 
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Xo(VSo)^[^(^)]„ = (^V^ - VSi • V)xo + S,xo, 

and therefore 

Xi = Xo / 7^[^(iv^ - VSi • V)xo + ^i]. (6.38) 

Near So = 0, we may use (6.36) to expand the first term inside the square bracket as a power 
series in A: 

— (^V' - VSi • V)xo = + 6_iA-i + bo + 

Xo 2 

biX + 62A2 + • • • , (6.39) 

where 6_2, 6-1, 60? bi, ■ ■ ■ are constants. Because of (6.37), in order that xi be analytic at 
q = 0, we must require 

Si = -bo; (6.40) 

otherwise, xi would have a term proportional to XolnA, which is not admissible. The function 
Xi is given by the integral (6.38) along the a-constant trajectory. In a similar way, we can 
determine £2, £3, ■ ■ • and obtain the solutions X2, Xs: • • •, in terms of quadrature along the 
classical (constant a) trajectory. 

7. PERTURBATION AROUND AN ATTRACTIVE COULOMB POTENTIAL 
Let He be the Hamiltonian for a Coulomb potential: 

^c = -^V^--, (7.1) 

where 

/ = Ze' (7.2) 
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with denoting the three-dimensional Laplacian, r the radius, Ze the nuclear charge and 
— e the electronic charge. Consider a problem in which there is an additional perturbation e\J; 
the corresponding Hamiltonian is 

H = H^ + eU, (7.3) 
where U is not singular at the origin. Let tpc and -ip be the ground-states of H^, and H; i.e., 

H,,!;, = E.^P, (7.4) 

and 

Hi) = Eifj. (7.5) 

7.1 Isotropic Case 

We first discuss the case that U{r) depends only on r. As we shall see, with modification 
our method can be adapted to derive ^ by quadratures along the radial trajectory. The solution 
to the Coulomb problem is well known: 

i)c = e-^"'- and = -^g\ (7.6) 

This suggests[10] that instead of (1.7)-(1.8), a different (7-power expansion is needed, one that 
should conform to the form ([Tfil ) of the Coulomb wave function. 
Write 

ij = e-^. (7.7) 

Expand 

S = g'So + Si+ g-^S2 + g-'S, + ■■■+ g-^^'^-^^Sn + ■■■ (7.8) 

and 
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E = g^Eo + g^E, + E2 + g'^E^ + ■■■ + g'^^^'-'^E^ + ■■■. (7.9) 

Since 

( [7. 51) becomes 

-hvSf + \v^S - — + eU = E. (7.10) 

^ Zi 'V 

Substituting (|7^ )-(|7^) into ( [TTUI) and equating the coefficients of g'^, g^, g^ , ■ ■ ■ , gf"^™, ■ ■ ■ 



on both sides, we obtain 



(VS'o)2 = -2Eo (7.11) 

V^o-VSi = Vso---i5;i (7.12) 
I r 

VSo ■ VS2 = -^(V5i)2 + iv^Si + eU- E,, (7.13) 

VSo • VS3 = -V5i ■ VS2 + ^V^S2 - Es, (7.14) 



etc.; for n > 1, 



n-l 1 1 

V5o ■ VS2n = - E ■ WS2n-m ' 7j(V5„)' + -V'^Sn-l - ^2n (7.15) 



2' ' 2 

m=l 



and for n > 1 



V^o ■ V^2n+l = - E ■ VS2„+l_^ + ^V^Ssn - ^Sn+l- (7.16) 

m=l ^ 

For the ground-state, these S'm(r) are all radial functions. Since U{r) is regular at r = 0, 
we can set 

U{0) = 0. (7.17) 

We adopt the same normalization condition (1.36), with 

7/^(0) = 1 and S(0)=0. (7.18) 
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From ( [7.111) , it follows that 



and (7.19) 



So = V^2E^ r. 



Substituting these expressions into (7.12), we find 



-dSi 1 



-2Eo— ^ = -(V-2^0 - 1) - E,. (7.20) 
or r * 

Because 5*1 should be regular at r = 0, y/—2Eo — 1=0 and therefore[10] 

Eo = -l and So = r (7.21) 
confirming ([7.61); ( p'.20[ ) then leads to 

and (7.22) 

^1 = -E, r. 



Since 



VV = -, (7.23) 



r 



we see that (7.13) becomes 



dSo 1 „o El s „ 

--El - — + eU{r) - E2. 



dr 2 
In order that S2{r) be regular at r = 0, 

El = 0, (7.24) 

and consequently, on account of ( [7.22[ ), 
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Si{r) = 



(7.25) 



and 

dr 

By using ( [7.23|) , we see that as r — > 



ef/(r) - E2. (7.26) 



V^Ss ^ (V252)o = -(^)o = -[eU{0) - E2]. (7.27) 

From (7.14), we see that in order to have S3 also regular at r = 0, (^)o must be 0; this 
yields on account of ( [7.17| ), 

E2 = eU{0) = (7.28) 

= eU{r) (7.29) 

and 



dSs _ e d . 2 



dr 2r^ dr 

As r 0, (7.14) implies 

883 .dSs. r e d . 2 



7r-2Tr{r'U) - E3. (7.30) 



dr dr 2r^ dr 

and, on account of (|7.23| ), 



r'f/)]o - ^3 (7.31) 



V^^3 - -(^)o. (7.32) 

r dr 

Hence by using ( [7.15| ) for 5*4 and in order that 5*4 should be regular at r = 0, 

(^) 

^ dr ' 

and therefore 



0, (7.33) 
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The same reasoning yields for all m > 1, at r = 0, 



dr 

From ( [7.34|) , we see that for 

eU{r) = er' 

E3 = 



|e if 1 = 1 



if / > 1. 



Example 



To carry the analysis further, let us consider the case / 

ef/(r) = er^. 

In this example, the above results give 

El = E2 = = 
Si =0, 82 = ^er' 

and 



Set n = 2 in ( [7.15| ), we have 



1 2 4 , Q TP 

— —-€ r + 3e — 



dr 2 

Because of (|7]3|), at r = 0, ^ = 0; hence, 

= 3e 

and 
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10 

In a similar way, we derive 



^4 = -^^eV. (7.40) 



^5 = 


0, 




^5 = -|6V, 




Eq — 


0, 


^6 


_ 43 2„3 I 1 ,3„7 
12 ~ 14 ' 




Ej = 


0, 




_ 43 2 2 I 13 3»,6 
4 12 ' 


(7.41) 


Es — 


129 2 

4^' 




etc. 





Putting together, for eU — er^, we find 



S = g^r+ -^r^ + ^r' - r^r^ 

JjLr^ + ±(_l^eV + — eV) 
8(?8 (^10^ 12 14 ^ 



+35(-T^''-' + T^^^''-'')+0(r[l) (7.42) 



and 



^ _^ + ^ _ i^e^ + 0(4). (7.43) 

In the usual perturbation series, to each order of the perturbation, the derivation of the 
perturbed Coulomb wave function requires summations over an infinite number of excited 
bound-states, plus the continuum. Here, to each order e™, the perturbed wave function can 
be obtained in closed form by quadratures along the radial trajectory. 

The perturbed energy can also be derived in an alternative way by using the integral form 
(4.92): 

oo 

E^—g'-r'-^^^ . (7.44) 



To first order in e, we can approximate the factor by e~^^^' in the integrals. For U = r^, 
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°° - 2 



^ e / e "^r^dr 

E = --g' + ^ + 0(a 



= -^^?^ + ^ + 0(6^). (7.45) 

To calculate E to the accuracy of 0{e^), we need the wave function to 0(e). From ([7.42|) , 

e-' = e-^\l - + ^) + ©(e^)]. (7.46) 



Substituting this expression into ( [7.44|) , we find the second term on its right-hand side to be: 

(7.47) 



6/e-^V[l-6(^ + ^)]dr 
" " 

a M 



confirming ( [7.43|) 



7.2 Stark Effect 

Consider the anisotropic example, in which the perturbation is 

eU = ercosa, (7.48) 

where a is the polar angle; i.e., 

= + + z'^ and z = rcosa. (7.49) 

Replacing ef/(r) by ercosa, we see that ( [7.1|) -( [7^28|) remain intact with Ei = E2 = Si = 0; 
(7.29) becomes 

— — = e r cosa, (7.5(J) 
ar 

and therefore 

5*2 = -e r^cosa, (7-51) 

and 
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V^5'2 = 2 e cosa. 



Thus, (7.14) gives 



— e cosa — E3, 



dr 

which yields, on account of (7.49), 

83^6 z- E^r 

and 



Since 5*4 satisfies (7.15) for n — 2, 

In order that 5*4 not have a Inr singularity, 

and therefore 

Ss^e z. 

In order that V'^S^ not have a 1/r singularity, 

E4 = 0. 

Thus, 



^4 = -—eV(l + 3 cos^ 
24 



Likewise, it is straightforward to derive 
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^5 


= 0, 


-^5 = ■ 


771 

-C/6 


9 2 


CI 

Se = 




= 0, 


67 = ■ 


Eg 


= 0, 




Eg 


= 0, 


5*9 = ■ 


Eio 


= 0, 


'S'lo = 


Ell 


= 0, 


^11 = 


El2 


355£ 
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7 

-— eV(l + cos2a), 
io 

= — eV^cosafl + cos^( 
16 

— e r cosald + cos a). 
48 ^ ^ 

— eVcosa —e'^r^i 

16 128 ^ 

53 3 99 4 4, 

— e rcosa e r (J 

8 512 ^ 



(7.59) 
(7.60) 
(7.61) 
(7.62) 
(7.63) 



761 

■— eV(l + 3cos2a)+0(6^), (7.64) 
31 31 

■^eV(l + cos^a)+0(e^), (7.65) 

(7.66) 



etc. 



Combining these results together, for the potential — ^ + ercosa, we find that, in powers 
of e, the wave function e"^ and the energy are given by 

er 1 e^r^ 7 1 

S^g^r + ^cosa(l + -/r) - ^[^^(1 + cos'a) + 24/^(1 + 3cos'a)] 

+ 0(e^) (7.67) 



and 



1 , 9 3555 



The — 1^ is known [11, 12]. Our method gives closed expressions for both the wave function 
and the energy to any finite order of e. 



68 



ACKNOWLEDGMENT 



One of us (T. D. Lee) wishes to thank T. K. Lee for discussions, and for a private communi- 
cation which outlines how our method can be extended to include the hydrogen atom. We also 
wish to give our thanks to RIKEN, Brookhaven National Laboratory and to the U.S. Depart- 
ment of Energy [DE-AC02-98CH 10886] for providing the facilities essential for the completion 
of this work. 



69 



REFERENCES 

[1] R. Friedberg, T. D. Lee, W. Q. Zhao, IL Nuovo Cimento A112(1999), 1195 

[2] A. M. Polyakov, Nucl. Phys. B120(1977), 429 

[3] E. Gildener and A. Patrascioiu, Phys. Rev. D16(1977)423 

[4] J. Zinn-Justin, Nucl, Phys. B192(1981), 125 

[5] V. A. Novikov, M. F. Shifman, A. I. Vainshtein and V. I. Zakharov, Usp. Fiz. Na 
Phys. Usp.) 136(1982), 553 

[6] A. A. Alejnikov and E. Shuryak, Yad. Fiz. (Sov. J. Nucl. Phys.) 46(1987), 122 

[7] E. Shuryak, Nucl. Phys, B302(1988), 621 

[8] C. F. Wohler and E. Shuryak, Phys. Lett. B333(1994), 467 

[9] S. V. Faleev and P. G. Silvestrov, Phys. Lett. A197(1995), 372 

[10] T. K. Lee, (Private communication) 

[11] H. M. Foley, 0. M. Sternheimer and D. Tycko, Phys. Rev. 93(1954), 734 
[12] A. Dalgarno and J. T. Lewis, Proc. Roy. Soc. A233(1955), 70 



70 



